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Abstract. Infinite energy solutions to the Navier-Stokes equations in R may 
be constructed by decomposing the initial data into a finite energy piece and 
an infinite energy piece, which arc then treated separately. We prove that 
the finite energy part of such solutions is bounded for all time and decays 
algebraically in time when the same can be said of heat energy starting from 
the same data. As a consequence, we describe the asymptotic behavior of the 
infinite energy solutions. Specifically, we consider the solutions of Gallagher 
and Planchon [l] as well as solutions constructed from a "radial energy decom- 
position". Our proof uses the Fourier Splitting technique of M. E. Schonbek. 



1. Introduction 

The purpose of this paper is to explore the large time energy decay in of 
solutions to the system 

(1.1) ut + u ■ Vu + Vp — Alt = —u ■ Vw — V ■ Vu, 

V • u = V • u = 0, 
u(0) = uo e L^{M?) 

where u is the velocity of an incompressible fluid, p is its pressure and is a specified 
external vector field satisfying 

(1.2) ||V%|U„(K2) < Ct-^-t + ^ 

for a = when 2 < 77 < oo and either a = or a = 1 when r/ = 00. Such a 
system arises naturally when considering infinite energy solutions of the Navier- 
Stokes equation, which includes the case of "rough" initial data in the plane. 
Recall that the Navier-Stokes equations for an incompressible viscous fluid are 

(1.3) wt + w • Vuj + Vp - Au) = 0, 

V • w = 0, w(0) = wq 

where w represents the velocity of the incompressible viscous fluid and p its pressure. 
The literature involving this equation is quite large and we mention quickly a few 
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relevant results. One of the first rigorous mathematical treatments of this system in 
the plane was the work of Leray [2] in which global existence of a unique solution 
corresponding to initial data in wq £ L^(IR^) was established. In questions of 
global existence and uniqueness are much more difficult and there are outstanding 
open problems even at the level of initial data. In however there has been 
much work dedicated to finding solutions with initial data in larger function spaces, 
for example see Gallagher and Planchon [I] , Cottet [3] , Giga, Miyakawa and Osada 
[1], Koch and Tataru [5], Germain [5], and the references therein. Particularly 
relevant to our discussion, in [1] and [6] the authors used a technique which involved 
separating the solution into a "rough" part and a finite energy part which satisfies 

Formally, if initial data wo is decomposed as wq = vq + mo with UQ e L^(R^) and 
if v{t) solves (jl.3p with initial data vq, then a solution of (|1.3p with data wq can 
be written as 'w{t) = u{t) + v{t) where u satisfies (jl.ip with initial data uq. The 
energy decay theorem we prove indicates that the energy of solutions to (jl.ip , that 
is |jM(t)||x,2(R2'), remains bounded and decays algebraically when the same can be 
said of the corresponding heat energy. In turn, this describes how w{t) approaches 
v(t) in the norm as time becomes large even though w and v need not be in 
individually. 

The main result in this article is the following Theorem: 

Theorem 1.1. Let u be a global solution to with v satisfying We 

assume: 

(i.) For each to > there is a constant Ct„ > such that for all T > to, 

(1.4) sup ||u(t)|U2(R2) <a„(i + r)i 

t,<t<T 

(ii.) For some 7 G [0,1] we have ||e^*uo||2 ^ C'(l + t)^^ where e^*uo denotes 
the solution to the heat equation with initial data uq. 

Then for every to > there exists a constant Ct^ such that 

\Ht)\\hiM^)<Ctoil + t)-'' 

for all t > tf). □ 

Remark 1.1. Assumption (i.) in the above theorem is the natural a priori energy 
estimate for (|l.ip . a formal proof is given in Subsection l2.1l Assumption (ii.) takes 
into account the natural decay rate for heat energy starting from ug. The Theorem 
states that if the heat energy starting from uq decays like (1 + 1)~'^ with 7 € [0, 1], 
then the solution u{t) of (jl.ip has the same energy decay rate. This is natural, as 
the heat equation is the linear part of (jl.l|) and we do not expect solutions to (|1.1|) 
to decay faster than this. On the other hand the "rough" terms (the nonlinear 
terms containing v) can "mix" the solution and slow the energy decay. □ 

Remark 1.2. It is known that the heat energy decay rate is determined by the 
behavior of ito near the origin in Fourier Space. For example, if mq € L^(R^) n 
L'^{R^), so that \uo{0\ < C near |^| = 0, then ||e^*wo||l < (1 + t)-^ and hence 
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ll'"(Olli^(R^) — ^(1 + 0^^- More detailed analysis may be found in Bjorland and 
Schonbek [7]. Although uq G i^(R^) implies ||e^*||2 — > as i — > oo, the heat 
energy may not decay at an algebraic rate (i.e. 7 = 0). This allows us to construct 
solutions to with arbitrarily slow decay by appropiately scaling the initial data 
and the external vector field, by using the same arguments as for the Navier-Stokes 
equations (for details on this case see Schonbek [8]). 

□ 

Remark 1.3. The proof of Theorem ll.il is based on the Fourier Splitting method of 
M. E. Schonbek [3], [S] introduced to study algebraic energy decay rates in parabolic 
equations. □ 

We now indicate how to use Theorem 11.11 to better understand the large time 
behavior of infinite energy solutions to 2D Navier-Stokes solutions. By an infinite 
energy solution we mean one belonging to one of the scale-invariant homogeneous 
Besov spaces Briq ^(K.^) which satisfy the chain of continuous embeddings 

(1.5) L2(R2) ^ b2A-1(m2) ^ ^2/r-i^j^2) ^ BMQ-^R^) C B^\^{^^) 

where 2 < r < f < oo and 2 < q < q < oo. 

Consider Navier-Stokes equations ()1.3p with initial data e S^.g' ^(R^), with 
r,q < oo. In these spaces, Gallagher and Planchon [1] proved existence of global 
solutions. To prove this result they decompose wq = vo + uq where uq S L^(R^) 
and vq € BMO^^IM.^) with small norm. Starting from this small vq they con- 
struct a solution v{t) of the Navier-Stokes equation using a fixed point argument 
which naturally satisfies ()1.2p for jy G [l,oo]. Next they consider (|1.1|) and find a 
solution u(t) using a fixed point theorem to obtain local existence then prove an a 
priori energy estimate to establish global existence. The energy bounds used imply 
||m(*)||l2(r2) < Ct^^^, which is exactly Assumption (i.) in Theorem 1 1.1[ though the 
authors leave open the question of finding better bounds on u. An interpolation 
argument is then used to show w{t) = v{t) + u(t) G Brlg'^i^"^) for all time. Us- 
ing similar methods, Germain [B] proved global existence of solutions for data in 
VMO'^iR'^) which is the closure of the Schwartz space in BMO^^iR"^). Moreover, 
he proved that under some mild conditions on r and q, Gallagher and Planchon's 
solutions with initial data in i?r,'q^~^(R^) are such that 

^hm ||u(i)||^v.-i(R2) =0. 

In this context we can use Theorem ll.ll to prove that the "finite energy part" of 
the infinite energy solution decays algebraically when the same can be said of the 
corresponding heat equation. This is the content of the following corollary: 

Corollary 1.2. Let wq G Br!q ^{M.'^), with r,q < oo. Consider wq = vq + uq, 
where uq G L^(R^) and vq G BM0~'^{'M?) with small norm. Let v{t) and w{t) 
be the solutions of ()1.3|) given in [1] with initial data vq and wq respectively. If 
lle*'^Mo||L2(K2) < C(l + ty^ for some 7 G [0, 1], then 
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\\wit)^v{t)U2^^2-^<C{l + t)-\ 

□ 

Remark 1.4- In particular, for any to > we have |lu||L°o([to,oo);L2(R2)) < oo which 
is stronger then the original energy estimate. □ 

More classically, Theorem 11.11 can be used to understand long time behavior of 
infinite energy solutions to the Navier-Stokes equations with finite local energy and 
integrable initial vorticity, that is wq € Lf^^{M.'^) and ljq = W x wq G _L^(K^). This 
initial data is a particular case of so-called vortex sheet initial data and it was used 
by DiPerna and Majda [10] to study approximate solution sequences for the Euler 
equation (see also [HI Sec. 3.1.2]). 

Remark 1.5. Initial data wq is of vortex sheet type if wq E Lf^^iM.'^) and loq = 
V X wq G A^(M^), where Ai{M.'^) is the space of nonnegative Radon measures. As 
for any ujq E A^(R^) there exists a unique wq G 

loq = curl Wo (see Corollary 4.4, Germain [6]), then wq is in one of the infinite 
energy spaces in (jl.Sp . □ 

Definition 1.3. An incompressible velocity field wq : ^ has a "radial energy 
decomposition" if there exists a smooth radially symmetric vorticity (Do(|a;|) such 
that 

\uo{x)\'^ dx < oo, 

where vq is defined from cZiq by the Biot-Savart law vq = K ljo, for K{x) = ^ |frp- 
the 2D Biot-Savart kernel. The radial energy decomposition, which is not unique, 
is possible on the whole plane since uq £ i^(R^) if and only if J^2 V x uq dx — 0. 

We restrict our attention to initial data with wq e Lf^^{M.'^) and ljq ~ W x wq E 
L^{R'^), because it can be split appropiately using the radial energy decomposition 
(see Lemma 3.2 in Majda and Bertozzi [H]). Moreover, some of the estimates we 
use when working with initial vorticity in L^(R^) need not be available in 7M(M^). 
Denote by ui{x,t) the solution to the heat equation with initial data ujq. As the 
initial data is radial, so is u}{x,t), and it is a solution to the vorticity formulation 
of Navier-Stokes equation 

(1.6) dtui + v -Vui = Aid, 

V = K * oj{t), 
Ldo{x) = w(a;, 0). 

With V ~ K * u) in hand we may then find the solution u{t) to p.ip starting from 
initial data uq using energy methods as outlined in [11] , thus obtaining the solution 
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'w{t) = v{t) + u{t) of the Navier-Stokes equation. In Subsection 12.21 we show how 
v{t) satisfies (|1.2|) . We have then the foUowing CoroUary: 

Corollary 1.4. Let w{t) be a solution of the Navier-Stokes equation with initial 
data wq G Lf^^{M.'^) such that ojq = V x ujq e L^(M^). Let wq = uq + vq be a radial 
energy decomposition with uq G L^(R^), wq a radial function, and vq — K * loq. If 
lle'^*uo|li2(R2) < C(l + for some 7 € [0,1], then 

lkW-«Wlli^(M^)<C(l + i)-^ 

where 

v{x,t) — -r—r^ / se^''a}o(s) ds 
fI Jo 

□ 



Remark 1.6. Using a far field calculation it can be shown that if V x has compact 
support then uq S LP(R^) for any p e (1,2] and 7 can be chosen to be any value in 
[0, 1). This is demonstrated in Subsection 12.21 □ 

Remark 1.1. For wq € ^^(M^), Gahay and Wayne [12], [13] have described the 
asymptotic behavior of solutions to the vorticity equation ()1.6p . In particular they 
prove 

X\^t^-\\\w{t) - (-^ ||l'!(R2)=0, 2<q<oo 

where a — ^0 dx and F(^) — Itt |fp (-'^ ~ ^ ■ Corollary concerns the 

borderline case q = 2, but we show how the solution approaches a radial solution 
instead of the Oseen vortex 

Note that the Oseen vortex is a solution to the Navier-Stokes equations (|1.3p with 
initial data wo('^) = ^ ifj?, which is not in Lf^^{M.'^), but is in i3r,oo~^(R^) because 
is a homogeneous distribution of degree —1 (see, Cannonc [Mj Lemma 3.3.2]). □ 

This articles is organized as follows. In the next Section we establish some basic 
properties of solutions to (jl.ip . including the a priori energy estimate. In Section 
[3] we use the Fourier Splitting Method to prove Theorem ll.il 

2. Preliminaries 

2.1. A Priori Energy Estimate. We now establish an a priori energy estimate 
for solutions of p.ip when v satisfies (|1.2p with a = and i] = 00. This estimate 
is known in the literature but we record it here for completeness since it is one of 
the assumptions for Theorem 11.11 It is straightforward to make it precise in the 
case of the radial energy decomposition mentioned in the Introduction (see [1] for 
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a rigorous argument in their setting). Formally, multiplying (jl.ip by u and then 
integrating by parts yields 

(2.7) + ll^"lli^(R^) <\<u-Vv,u>\ 

where we have introduced the notation < w- Vd, u >= / u-'VviUi dx. Fix to > 0. 
After integrating by parts and using Holder's inequality, then (|1.2p with a = and 
?/ = oo and then Cauchy's inequality, we have for any t > to, 

I < u{t) ■ Vw(t), > I = I < u{t) ■ Wu{t),v{t) > I 

In the above line the constant may depend on to- Combining this estimate with 
p.7p and then integrating from to to t yields 

\\u{t)\\l2^j,2^ < C f h(s)||i2(R2)(l+s)-lds+||^i(to)||i2(R2). 
Jto 

From here a Gronwall inequality gives 

\\umhiw2)<C\\uito)\\h(t2^il + t) 
which is (|1.4p of assumption (i.) in Theorem ll.il 

2.2. Properties of solutions with the radial energy decomposition. In this 
subsection we consider the Navier-Stokes equation with initial data wq G L^q^,(K^) 
such that ujq = W X wq ^ L^(R^). As in the Introduction, consider the radial energy 
decomposition wq ~ uq + vq where uq G L^(]R^) and vq is the velocity of the radial 
vorticity ujq. We first we prove our claim that v = K * e^*cZio satisfies the estimate 
(|1.2p with a = and i] = oo. As wo G i^(M^) we have by direct calculation 

||e^*wo||LP(M^) < Ct-(i-^), l<p<oo. 

To find the estimate on v{t), the corresponding solution to the Navier-Stokes equa- 
tions, we recall the following estimate on the Biot-Savart Kernel. 

Lemma 2.1. Let ljq G LP{R^) n L9(M2) forl<p<2<q<oo and let < a < I 
be such that -i = — + ^—^ . For v = K * cjo we have 

z p q 

||i'(0IIl°°(E2) < C\\ujQ\\1p^jg.2-^\\LJo\\]~^^2y 

Proof. See [H Lemma 2.1]. □ 

Combining the previous lemma with the above bound on e^*wo we establish 
(|1.2[) . As mentioned in the Introduction, if we further assume that ujq = 'V x uq has 
compact support Bn we can use a far field calculation to demonstrate ||e'^*tio||2 < 
C(l + t)~'' for every 7 G [0, 1). Indeed, ii y < R and x > 4i? then the following 
geometric series converges: 

1 _ 1 ^f\y\^ 2x-y\' 
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Using /jja V X uo rfic = we find that for large x 

"o(a;) = ^ f--r^ / y-^uJo{y)dy - -^x ■ [ yoJo{y) dy + 0{\x\~^) 

which implies that uq G L^(M^) for every p e (1,2]. For q such that | + ^ = f we 
bound 

||e^*-tio||2 < ||'S(i)||L<!(R2)||wo||L<j(R2) 

where <i>(t) is the 2Z? heat kernel. As ||$(t)||ig(R2) < Cf^"^ it must be that 
||e^*uo||2 < Cti^p. Since p e (1,2] we have |ie^*Mo|li < C{1 + t)-^ for every 
7 6 [0,1). 

3. Decay 

In this section we prove Theorem 1 1 . 1 1 using the Fourier Splitting Method of M. 
E. Schonbck. In our proof we also incorporate a Gronwall-type trick used by Zhang 
|15j . Here we proceed formally but note the argument can be made rigorous using 
an approximating sequence of solutions. This would be argued similar to the proof 
of the energy inequality (|1.4p in [T] or similar to [5] in the more classical radial 
energy decomposition case. We start with frequency bounds. Applying Duhamel's 
formula in Fourier space and a simple integral inequality to (|l.ip yields 

(3.8) \u\ < e-l«l'*|wo| + ^ e-l«l'(*-^)|C| + l^iTSl + + ds. 

Taking divergence of (|l.ip and then using the symmetry of the tensor we find that 
\p\ < 2\v (S) u\ + \u (E) u\, so we obtain the bound 

|u| < e-l«''*|uo| + 2 e-l«l'(*-")|e| + I^Tm^l) ds. 

Using now Holder's inequality with the estimate (|1.2p (r] ~ oo and a = 1) gives 

\<u-S/v,u>\<Ct-'\\u{t)\\l,^^,^ 
so that after multiplying the PDF by u and integrating by parts we have 

We fix to > and now consider the inequality for < > to > so that < 
(l+V^)(l+t)-i and 

Ct-l||u|li2(R2) <Co(t+l)-lh||i2(R2), 

where Co contains the term (1 + tp ^). We now apply a Fourier Splitting argument 
around a ball with radius r{t) > 0, where r{t) is to be determined later. After 
observing that 
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we find that 

(3.9) J4kl!i2(E^) + (^'-Co(i + l)-^)hlli.(«.)<r^ / \u{s)\'d^ 



B(r) 



for t > to. 

In the case where (|1.2p does not hold for = oo, a = 1 we can instead use (|1.2p 
with rj = oo,a = as mentioned in the Introduction. After integration by parts 
and using Cauchy's inequality we obtain the bound 

\<U- S/V,U > I I < W • WU,V > I < ||u||i2(R2)||VM||i2(R2)||i;||ioo(R2) 



1, 



<C(l + i)-^||«(i)||i2(R2) + -IIV7.il 



2 

LH 



Considering again a fixed to > we again arrive at p.9p but with different constants 
which will have no impact on the following arguments. Thus, we can say that 
Theorem 1 1 . 1 1 holds for these two estimates on v, which is what we use as hypotheses 
in our Corollaries 11.21 and 11.41 

Now we estimate the right hand side of p.Q^ : 

\uis)\^di< I e-2|«l^*|Ko|2de 

B{r) JB{r) 

+ [ ([ e-l«l'(*-^)|e|(|7T^| + |^I^|)ds) 

JBir) \Jo J 

■.= I{t) + B. 

We need to break B into two pieces, one with |m (g) u\ and the other with \u ^ v\. 
This is done with Minkowski's inequality then the triangle inequality by 



B < 



B(r) 



f e-l«l'(*-«)(|,7^| + |S:®li|)<is^ d^ 



<r^\l U e-2|«l'(*-^)(|,^| + |ir^|)2d^j ds 



e-2|«l'(t-«)(|,}^|)2d^] ds+ f i I e-2|?l'(*~^)(|u(^u|)"de 

\JB{r) J Ja \JB(r) y 

Using Holder's inequality, then the decay assumption on v (here rj ^ oo) 



e 



-2i?i'(t--)n^i^|)2d^ < / e-^'i«i'(*-^) I II^^®"IIl.(r2) 



B(r) J \-'B(r) 

< C(t - s)"p ||-u(s)||i2(R2)||i;(s)||ip(R2) 

< C(t- s)"p||M(s)||i2(R2)S-(5-^) 



DECAY OF INFINITE ENERGY SOLUTIONS TO NAVIER-STOKES EQUATIONS 9 
2 9 ^ p- 



where in the above sequence 5 = ^ + ^- Also 



e 



B{t) 



<C|r|||«(,s)l|i.(R.) 



so that 



Then p.9p becomes 



+ Cr^( f \\u{s)\\l,,^..ds^ 



Choose r {t) ~ jt+i) ^11^ multiply everything by 2 (t + 1) to find 









) i2(E2) ds^ 





By assumption we have 

(3.10) I{t) <C{l+t)-'^ 

for some 7 € [0,1]. The next step is to integrate from Iq to p and divide by 
(1 + pY~'^ ; which leads to 



{l + pV\\u{p)\\l.^^. 



(1 




(H 






c 


(n 






c 


(H 





2 



c 



L2(R2 



) + (rT7j^ Jo + + 
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The main goal now is to set it up as a Gronwall inequality for g{p) = (l+p)'''||u(/9)||2. 
For the Ai term we have 

* 1 11 

{t - s)"P ||w(s)||l2(r2).s"(2-p) ds 

< (^j\t~s)-h-'''-i\i+s)-usy (^j\i+s)\\u{s)\\h^j,,)dsy 

<c(^l\l + s)\\u{s)\\l,^^,^dsy . 
Here we used J^{t — s)~ps~^2-?) ds < C for alH > when p > 2. Then, 



Jo 

In moving to the last line we used the fact that 7 < 1- The A2 term is similar, as 



< 



< 



c 








c 


(1- 






c 


(1- 


f p)2-7 



2 



Now, by the assumed bound (|1.4p . ||m(s)|||2(r2) < C(l + s) so that 



Jo 

The term (jqi^p^ll'^l^o)!!! bounded by some constant. Using the assumption on 
I{t) 

c r\ w ^ , c 



(l+p)2"7 

Putting everything together we have 

9{p)<C + C r g{s)ds 
Jo 

9{p) = {l + pV\\u{p)\\l 
so Gronwall's inequality implies g{p) < C or 

\\u{t)\\l<C{l + t)-\ 
This is exactly the conclusion in Theorem ll.il 
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